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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2411 

METHOD OF DETERMINING INITIAL TANGENTS OF CONTOURS OF 
FLOW VARIABLES BEHIND A CURVED, AXIALLY 
SYMMETRIC SHOCK WAVE 
By George P. Wood and Paul B. Gooderum 


SUMMARY 


A general method for calculating the initial tangents of contours 
of constant density, pressure, and Mach number behind a curved, axially 
symmetric shock wave has been derived. Initial tangents of density con- 
tours are obtained that fair -well into contours determined from inter- 
ferograms for the flow about a sphere at Mach number 1.62. Streamlines 
and contours of constant Mach number throughout this flow field have 
also been deduced from contours of constant density obtained by inter- 
ferometry. The flow field obtained from the shock wave by the method 
of characteristics is compared with that obtained by interferometry. 


INTRODUCTION 


Various phases of the problem of flows behind curved shock waves 
have been treated theoretically by many investigators. Crocco, in a 
pioneer paper (reference 1), investigated two-dimensional flow behind 
a curved attached shock wave. Significant contributions were made later 
by Maccoll, Drougge, and Guderley. Maccoll (reference 2) calculated 
the flow around various bodies with detached shock waves at near- sonic 
velocities. Drougge (reference 3) calculated the pressure distribution 
on conical tips with detached shock waves. Guderley (reference 4) 
investigated the transition between flow with a detached shock wave and 
flow with an attached shock wave. Thomas, in a series of papers of 
which reference 5 is representative, considered the curvature of attached 
shock waves and the pressure distribution on bodies behind attached shock 
waves. Lin and Rubinov (reference 6) obtained general relations for flow 
behind curved shock waves and developed a method for calculating the sub- 
sonic flow between the shock wave and the body. This method was applied 
by Dugundji (reference 7) to calculate the pressure distribution along 
the axis between a sphere and a shock wave. Busemann (reference 8) 
reviewed analytical methods for the treatment of two-dimensional flows 
with detached shock waves and discussed the main features of mixed 
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subsonic and supersonic flows. Ferri (reference 9) presented a method 
for obtaining the pressure drag in two-dimensional or axially symmetric 
flow behind a detached shock wave when the location and shape of the 
shock wave are known. Moeckel (reference 10) gave an approximate method 
for predicting the form and the location of detached shock waves. 

Experimental investigations of the flow behind curved shock waves 
include the investigation by Ladenburg, Van Voorhis, and Winckler 
(reference 11) who determined the density contours around various axi- 
ally symmetric bodies at Mach number 1.7* The present authors in refer- 
ence 12 presented data on the shape and location of detached shock waves 
on cones and spheres over a range of Mach numbers from 1.17 to 1.8l. In 
reference 13 interferograms of the flow around a sphere at Mach num- 
bers 1.30 and 1.62 were evaluated; the contours of constant density 
ratio in the flow field around the sphere were determined. However, the 
initial tangents of the density contours could not be obtained accu- 
rately from the interferogram by the method of analyzing interferograms 
that was used in reference 13 . 

The purpose of the present paper is to present a method for calcu- 
lating the initial tangents of the contours of flow variables immedi- 
ately behind a curved, axially symmetric shock wave when the location 
and shape of the shock wave are known. This method was used in the 
present paper to obtain the initial tangents of the density contours in 
the flow field around a sphere at Mach number 1.62. Streamlines and 
contours of constant Mach number in the flow behind the shock wave have 
also been obtained from the density contours of reference 13 . The flow 
field in part of the supersonic region behind the shock wave has been 
computed by the method of characteristics (reference 14) and the results 
compared with those obtained from the interferogram. 


SYMBOLS 


a velocity of sound 

A cross-sectional area 

b coordinate in plane of binormal 

Cp specific heat at constant pressure 

c v specific heat at constant volume 


D 


diameter of sphere 
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F 


log 

M 

n 

P 

R 


Rw 

S 

AS 

t 

T 

T b 


T n 

V 

x 

y 

7 

5 

e 

% 


distance to focus of streamline directions due to curvature 
of shock wave 

logarithm to base e 

Mach number 

coordinate normal to direction of streamline 
pressure 

radius of curvature of streamline 
radius of curvature of shock wave 
entropy 

change in entropy - Sq^ 

coordinate tangential to direction of streamline 

total taper distance of element of stream tube 

taper distance of element of stream tube in plane of 
binormal 

taper distance of element of stream tube in plane of 
normal 

velocity 

coordinate in direction of axis of symmetry 
coordinate normal to axis of symmetry 
ratio of specific heats ^Cp/c y ^ 

angle of deviation of flow across shock wave 
angle of shock wave 

angle between direction of x-axis and initial tangent of 
contour of constant Mach number 

angle between direction of x-axis and initial tangent of 
contour of constant pressure 
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angle between direction of x-axis and initial tangent of 
contour of constant density 


IT 


ratio of stagnation pressures 



P 


density 


a 


coordinate along shock wave 


Subscripts : 


s 


stagnation 


0 


free stream 


SHOCK WAVE 


In order to apply the method of calculating the initial tangents, 
that is developed in the present paper, the location of the shock wave 
and the variation along the shock wave of both the slope and the radius 
of curvature along the shock wave must be known. When the location and 
the slope variation of the shock wave are known, the streamlines can be 
obtained and the method of characteristics can be used to determine the 
flow field. The present section, therefore, is a discussion of the shock 
wave, its first and second derivatives, and its radius of curvature. 

Figure 1, which is taken from figure 4(b) of reference 13, is the 
interferogram from which the coordinates of the bow shock wave and the 
other experimental results are obtained. This interferogram shows the 
flow about a sphere at a free-stream Mach number of 1.62. The free- 
stream stagnation temperature was 533*1° F absolute, and the free-stream 
density was 0.003501 slug per cubic foot. The rim shock wave from the 
edge of the nozzle is shown in the upper left-hand corner of the inter- 
ferogram. This shock wave is very weak, as the fringe displacement 
across the shock wave is very small and the shock wave lies within 1° 
of the Mach angle. Because of the rim shock wave, however, and the near- 
ness of the edge of the jet, the portion of the bow wave of the sphere 
between the rim shock wave and the edge of the jet was considered to be 
not sufficiently accurate and was not used (except for calculating one 
of the initial tangents; namely, the contour of constant density ratio 
of 1.5). The location of the portion of the bow shock wave between the 
rim shock wave and a position near the axis of symmetry was carefully 
measured and is shown in figure 2, where y is the normal distance from 
the axis of symmetry, x is the distance along the axis from the inter- 
section of the shock wave and the axis, and D is the diameter of the 
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sphere. Because of the difficulty of making very accurate measurements 
near the axis where the slope of the shock wave is very large, the part 
of the shock wave near the axis has been shown as a dashed line. 

The variation of the slope of the shock wave is given in figure 3* 
This curve was obtained by first reading y at values of x that 
increased by steps of 0.1, and then determining Ay, the increment in y 
for the given increment Ax in x. Next Ax/Ay was -plotted against y. 
The value of x was computed by finding the area under each small por- 
tion of the curve. If the value of x so obtained checked closely with 
the measured value of x, within approximately one part in a thousand, 
then the curve of Ax/Ay was taken to be correct. If the check was not 
close, the curve was refaired until a close check was obtained. Then 
points were read from the curve of Ax/Ay and were used for plotting 
the curve of dy/dx shown in figure 3* 

The variation of d 2 y/dx 2 was obtained in a similar manner. The 
Ax Av 

quantity — ; was plotted against — and the curve was refaired 

A(Ay/Ax) Ax 

until integration gave the original value of x within approximately 

I d^y \ x 

three parts in a thousand. The variation of -I — -ID with =- is 
shown in figure 4. yrj D 

From the variation of the first and the second derivatives, the 
variation of the radius of curvature of the shock wave was calculated. 

The radius of curvature obtained by this method, however, is not very 
accurate. Errors in the location and in the slope of the shock wave are 
magnified in the second derivative and, therefore, in the radius of curva- 
ture. In order to reduce the uncertainty in the radius, the shock wave 
on the other side of the axis (most of which is not shown in fig. 1, but 
which is included in the original interferogram from which fig. 1 was 
taken) was measured and its derivative and radius were calculated and 
plotted. The radius agreed fairly well with that obtained for the other 
side of the axis. Furthermore, with a mechanical device and a very large 
reproduction of the interferogram, the normals to the shock wave, the 
envelope of the normals, and then the variation of the radius of curva- 
ture along the shock wave were obtained. These results agreed rather 
well with the calculated values. A curve showing the variation of the 
radius of curvature along the shock wave that results from fairing 
through the measured points and the calculated points on both sides of 
the axis is given in figure 5- 
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RESULTS AND DISCUSSION 

Initial Tangents of Contours of Flow Variables 


In reference 13 the interferogram of figure 1 was analyzed to 
obtain the contours of constant density behind the shock wave. These 
contours are shown in figure 6, which is taken from figure 5(b) of 
reference 13. As was explained in reference 13, the analysis did not 
give accurately the contours in the region immediately behind the shock 
wave. For that reason there was a gap in figure 6 between the beginning 
of each contour and the known theoretical point where the shock wave and 
the contour should intersect. The purpose of the present section is to 
fill in these gaps by giving the theoretical values of the tangents of 
the contours at the shock wave. 

The analysis for finding the initial tangents of the flow vari- 
ables is given in detail in appendix A; however, the method is 
described briefly in the present section. The pressure immediately 
behind a curved shock wave varies along the shock wave. The rate of 
variation is obtained by differentiating along the shock wave the well- 
known equation for the pressure ratio across an oblique shock wave. The 
general rate of pressure change in any direction at any place in the 
flow field - including the one along a shock wave - may be represented 
by the two basic pressure derivatives, namely, in the direction of the 
tangent to the stream tube and in the direction of the normal to the 
stream tube. Therefore, the rate of pressure change along the shock 
wave is (appendix A) 

S-iHS-E = cos( e - 5) Ll2S_E + 6ln(6 . 6) a log p (1) 

dc dt dn 

Equation (l), with known values of e, 8, and — is one rela- 

d(T 

tion between the two pressure derivatives. Another relation is also 
needed. 

The tangential rate of change of pressure results from the tangen- 
tial rate of change of velocity that is necessary to meet the require- 
ment of continuity when there is a change in the cross-sectional area 
of the flow, and is given by Newton's second law of motion as 
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or 


d log p _ 7 M 2 d log A / 

~^ = 1 

The other basic rate of change is normal to the flow and balances the 
centrifugal force 

dp = pV 2 
dn " " R 


or 


d log p _ 7 M 2 

dn R 


( 3 ) 


An equation that relates the two pressure derivatives can be 
obtained from equations ( 2 ) and ( 3 ) by finding the relation between 

— A and R, as follows: At every point on the shock wave, the 

at 

local slope of the shock wave determines the initial streamline direc- 
tion behind the shock wave. Inasmuch as there is always a point at the 
shock wave where the streamlines have no initial curvature, that point 
may be used for orientation, and the initial tangents of the streamlines 
may be considered rather than the streamlines themselves. Because the 
shock wave is curved, successive initial tangents of the streamlines 
behind the shock wave are not parallel and, therefore, if extended, con- 
verge to a focus (or in larger extent to an envelope) in the focal dis- 
tance F (fig. 7)* If the streamlines were straight, the focal dis- 
tance F would determine the rate of change of cross-sectional area of 
the stream tubes and the rate of pressure change along the stream tube. 
Contours of constant pressure would be concentric circles around the 
focus (or the evolutes of the envelope). Since the pressure at the 
shock wave is known and since the entropy is known and is constant along 
a stream tube, there are two independent ways of calculating the rate of 
pressure change along the stream tube, one from the given pressure con- 
tours and the other from continuity inside the stream tube. These two 
methods, however, will agree for the tentatively straight streamlines 
only when the curvature of the shock wave and the distance from the axis 
of symmetry happen to have, fortuitously, the necessary values for agree- 
ment. In all other cases, a finite initial curvature of the streamlines 
has to modify the rate of change of area and to shift the contours of 
constant pressure in such a way that the two methods give the same 
result. In figure 8 the initial curvature of a streamline is indicated. 
In appendix A the details are given of finding the relation between the 
initial radius of curvature R of the streamline and the rate of change 
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of area of the stream tube. This relation is then used with equa- 
tions (2) and ( 3 ) to obtain the following equation between the two pres- 
sure derivatives: 



+ c sc(e - 6) — (4) 


da 


Equations (l) and (4) form a system of two equations in two 
unknowns and can be solved for the unknowns, the pressure derivatives. 
Then the effect of entropy variation along the shock wave can be taken 
into account to obtain expressions for the density derivatives and the 
Mach number (or temperature) derivatives. 

When the two derivatives of a flow variable are known at a point 
on the shock wave in perpendicular directions, they can be added vecto- 
rially to find the direction of the maximum rate of change of the vari- 
able. Normal to that direction, then, is the direction of zero rate of 
change, which defines the tangent of the contour of constant value for 
the variable under consideration. 

The final results are given by the three following equations, in 
which 9 is the angle between the direction of the x-axis and the ini- 
tial tangent of the contour of constant value of the variable indicated 
by the subscript on 9 (the angle 9 is considered positive when meas- 
ured counterclockwise from the direction of the x-axis): 


(tan ©p - 5) 


sin(e - 5) sin 5 + cot(e _ 5 ) 
_ ?M 2 


8 log p dS 
+ de 


tan(e - 6) 



8 log p 85 


8e 


( 5 ) 


to 


tan(0p- 5) 
tan(€ - 6) 


— sin(6-5)sin 5 + — = ^ cot (e - 5) ^ B + 
y 7JT ° € 

5g sln(6-5)sln 6 + »^U tan(€-8) 5 ^ _ r-j, 1 - (m 2 - l)tan 2 <6 -5) 8 l o g II 

y 7M 2 & 7 M 2 tan(€ - 5) 8* 

( 6 ) 


tan ( e M ~ S ) 

tan(e - 6) 


R. 


w 


sin(e - 6) sin 5 + — — cot (e - 5) 


yVf 


— sin(e - &)sin 5+ — ~ tan(e 

y 7 i £ 


, 5 ) d log p t 5& [ 1 - 
8c 8e 


log p 85 
8e 8e 

(m 2 - l) tan 2 (e -6) 8 log IT 
7W^tan(e-6) de 


(T) 


The points on the shock wave where the contours of constant density should intersect the 
shock wave are known from the shock-wave angle. The directions of the tangents of the density 
contours have been calculated for these points by equation (6). The results are shown in the 
following table, which gives the value of the density ratio for a contour, the point on the 
shock wave where the contour should meet the shock wave, and the calculated values of 0p at 
that point : 


Density 

ratio, 

p/p 0 

x/D 

e P 

(deg) 

2.0 

0.046 

64 

1.9 

.097 

39 

1.8 

.168 

ll 

1.7 

• 309 

-21 

1.6 

.508 

-70 

1.5 

.823 

-95 
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The tangents are also indicated in figure 6 by the short, straight lines 
that have been drawn through the shock wave. The initial tangents given 
by the analysis fit well between the experimentally determined contours 
and the known points where the contours should intersect the shock wave. 


Streamlines and Contours of Constant Mach Number Determined 

from Density Contours 

The only known quantity in the flow field behind the shock wave, 
except at the shock wave and at the surface of the sphere; is the den- 
given in figure 6. The values of the other state variables of 
the gas are not immediately known because of the variation along the 
shock wave of the entropy behind the shock wave. The streamlines, as 
the contours of constant entropy, can be found, however, by a graphi- 
cal method that is based on satisfying continuity along a stream tube. 
The local cross-sectional area of a stream tube can be expressed, if 
its entropy is given, in terms of known quantities in the free stream 
and the local values of density and stagnation density. The first 
stream tube is constructed by taking the axis and the surface of the 
sphere as a streamline. Local values of density, as given by figure 6, 
and the average value of stagnation density for the stream tube are used 
to calculate local values of cross-sectional area. Then, adjacent 
stream tubes are added in sequence. The method is described in detail 
in appendix B. 

From the steamlines and the density contours, the other flow vari- 
ables can be found. In the present case, the contours of constant Mach 
number were found. The results are shown in figure 9> where the density 
contours taken from figure 6 are shown by the heavy solid lines, the 
streamlines, by the light solid lines, and the calculated Mach number 
contours, by the dashed lines. 


Streamlines and Contours of Constant Density and Mach Number 
Determined by Method of Characteristics 

Contours of the flow variables in part of the flow field behind 
the shock wave were calculated by the method of characteristics. The 
part of the flow field for which results can be obtained is quite 
limited in extent. The flow field is limited on one side by the sub- 
sonic flow in front of the sphere and on the other side by the fact that 
the location and shape of the bow shock wave are not known accurately 
enough for use beyond the nozzle-rim shock wave. The method used was 
that given by Ferri (reference lA) . The only data needed for the calcu- 
lations were the free— stream conditions, the coordinates of the shock 
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wave (fig. 2), and the variation of the slope of the shock wave (fig. 3). 
The calculations were made with the Bell Telephone Laboratories X-66744 
relay computer at the Langley Laboratory. An attempt to obtain good 
accuracy in the results was made by using fairly small increments 
(Ay on the shock wave approximately equal to 0.04D) to give a network 
of characteristics that had a small mesh and by reiterating each step 
of the process. The results are shown in figure 10. Contours of con- 
stant ratio of density to free-steam density are shown in figure 10(a), 
where the contours obtained by interferometry are also reproduced from 
figure 6. Contours of constant Mach number are shown in figure 10(b), 
which also shows the Mach number contours obtained by satisfying conti- 
nuity, transferred from figure 9« Streamlines, or contours of constant 
entropy, obtained by the method of characteristics are shown in fig- 
ure 10(c), together with those taken from figure 9* 

The contours obtained by the method of characteristics and the con- 
tours obtained from the interferogram do not agree very well. Actually, 
however, the agreement is perhaps as good as might be expected. In 
deriving an axially symmetric flow field from experimental measurements, 
results that have an extremely high degree of accuracy are difficult to 
obtain. In the present case, the results obtained by the method of 
characteristics have some small percentage of uncertainty because of the 
finite size of the mesh of the characteristics network. The analysis of 
an interferogram of axially symmetric flow must also be made very care- 
fully in order to obtain results in which the uncertainty is very small. 
In reference 13, from which the interferometric results were taken, the 
estimate is made that the uncertainty in the location of any contour 
does not exceed plus or minus approximately one-half the distance 
between the adjacent contours. This amount of uncertainty, which is 
only about 3 percent, is just about equal to the discrepancy between the 
results obtained by the two methods. 

The fact that the streamlines are straight in the only region of the 
flow in which the method of characteristics could be applied is entirely 
fortuitous in the sense that the curvature of the shock wave and the 
distance to the axis for this region are such that the streamlines are 
straight. 


CONCLUDING REMARKS 


A method of calculating the initial tangents of contours of con- 
stant density, pressure, and Mach number at a curved symmetrical shock 
wave has been derived. The method gives tangents of density contours 
that fair well between the contours obtained by interferometry and the 
known points of intersection on the shock wave. 
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The streamlines of the flow about a sphere at Mach number 1.62 have 
been deduced from experimentally determined density contours by satis- 
fying the requirements of continuity along a stream tube. 

Part of the flow field around the sphere was also obtained from a 
knowledge of the shock wave by applying the method of characteristics. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , April 12, 1951 
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APPENDIX A 


CALCULATION OF INITIAL TANGENTS 


OF CONTOUES OF FLOW VARIABLES 


The problem is to find the tangents of the contours of constant 
pressure, density, and Mach number immediately behind an axially sym- 
metric curved shock wave. The solution given here was obtained by using 
ideas and suggestions of Dr. A. Busemann of the Langley Laboratory. 

A curved, axially symmetric shock wave in a uniform, steady, super- 
sonic free steam is considered. The location of such a shock wave can, 
because of axial symmetry, be completely given by two variables, and 
the shock wave can be drawn as a curve in a meridian plane. The local 
slope of the shock wave with respect to the free-stream direction of 
flow, together with the state of the gas in the free stream, completely 
determines the initial state of the gas and the initial tangent to the 
direction of the streamline at any and every point on the downstream 
side of the shock wave. In particular, the pressure is given by the 
equation 


Because the shock wave is curved, the local slope, the flow variables, 
and the initial tangent to the streamlines vary along the shock wave. 
In particular, the rate of change of pressure along the shock wave is 
obtained from equation (Al) as 


Another way of viewing the pressure change along the shock wave is 
the following: The rate of pressure change in any direction at any 
point in a flow field - for example, along a shock wave - is the sum of 
the components in that direction of the two basic rates of change. One 
of these basic rates of pressure change lies in the direction tangent 
to the direction of flow. This tangential rate of change results from 
any tangential rate of change of velocity that may be necessary to meet 
the requirements of continuity when there is a change in the cross- 
sectional area of the flow. It is given by Newton *s second law of 
motion as 



(Al) 


d log p _ de 2yMo g sln e cos 6 


(A2) 


yMo^sin^e - ^ 



Ik 
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By using the relations 


m 2 . v 2 |e 

dp 


and 



the equation can be written as 


5 log p _ y M 2 S log A 

St l _ m 2 St 


(A3) 


The other basic rate of pressure change is normal to the flow 
direction, balances centrifugal force, and is 

Sp _ pV 2 
Sn " R 

or 


S log p _ yM 2 
Sn B 


At any point on the shock wave, the rate of change of pressure along the 
downstream side of the shock wave is the sum of the components along the 
shock wave of Sp/St and Sp/Sn, that is, 


j ° g P = cos(€ - 5) — g - p + sin(e - 6) 

St Sn 


(A5) 


Another relation between - and — p i s needed. This rela- 

St dn 

tion can be obtained by finding first a geometrical relation between 
— - and R, the unknown factors in equations (A3) and (A4). Con- 


sider first a section of the flow as seen in a meridian plane (fig. 7) 
Because the shock wave is curved, the initial tangents to the stream- 
lines behind the shock wave change along the shock wave. The initial 
tangents, therefore, if extended, converge to an envelope or focus, in 
the focal distance F. In the general case, the stream tube is not 
identical with the extended initial tangents that converge in the dis- 
tance F, since the stream tube is curved in the meridian plane. If, 
however, the streamlines were straight, the distance F to the focus 
would be one of the factors that determine the rate of change of area 
of the stream tube, as this distance would define the taper of the 
stream-tube area in the meridian plane. The taper of a stream tube as 
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seen in a meridian plane is, however, not the only taper. Consider an 
element of an annular axially symmetric stream tube as shown in fig- 
ure 11. The cross-sectional area of the element is 

dA = dn db 

The element tapers not only in the meridian plane, which is the plane 
of the normal, but also in the plane of the binormal. In other words, 
both dn and db vary along the stream tube. The quantity db is 
proportional to y as shown in figure 11. If the cross section under 
consideration were moved away from the shock wave, y and db would 
become larger. If the cross section were moved in the opposite direc- 
tion, db would become smaller at a rate that would make it vanish at 
the axis of symmetry. The side of the element in the plane of the 
binormal therefore would taper to an edge or focus at the axis in the 

distance -Tv, = — — — 

D sin 6 

The total taper of the element of stream tube is made up of the 
components in the plane of the normal and in the plane of the binormal. 
The total taper of the element could then be considered as having an 
envelope distance T, and would be given by 

T F T b * A dt 

(where straight streamlines in the meridian plane are assumed). 

If the local Mach number in the element of stream tube were unity, 
the cross-sectional area of the element would be constant, dA/dt would 
be zero, and the normal and the binormal tapers would be opposite and 
equal ^F = At a Mach number of 1 behind the shock wave, however, 

d6/dcr is almost equal to zero for a perfect gas and the distance F is 
not equal to -T^. In fact, F is greater than -T^. In order, there- 
fore, for the taper in the plane of the normal to have the correct value 
of 


-T* = £ - 

^ sin 6 

the streamline must have an initial concave curvature at the shock wave, 
as shown in figure 12. The initial tangent on one side of the stream 
tube and the tangent at the opposite point of the stream tube, which is 
at a distance dt from the shock wave, are extended and intersect at 
the taper distance T n . For a Mach number of 1 the requirement is, 
therefore, that 


T n = -Tb 


1 6 
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and the streamlines must have a curvature to meet this requirement. In 
the general case, also, of Mach number other than unity, the envelope 
distance given by 


F = 


d n 
d6 


at the shock wave is not the correct taper in the plane of the normal, 
but is either too large or too small, and the streamline must have an 
initial concave or convex curvature. The component of taper in the 
meridian plane is defined by a different pair of tangents from the ini- 
tial tangents. The taper is always T n , as shown, but in the general 
case T n is not equal to -T^. 


The relation between the radius of curvature R and the taper 
distance T n can be obtained from figure 12, as follows: 6 varies 

both along cr and along t. The angle subtended by dt is — dt. 

dt 

The angle subtended by dn is - — dcr + — dt. The two triangles 

da dt 

with — dt as vertex angle are similar, therefore 


dt 

R 


dn + T n da 
n da 


T 


n 



d6 


or 

^ = tan(e - 5) i + sec ( e - 6) JJ 

This relation can now be used to obtain the relation between l/R and 
d log A 

St ' 


Equation (A 6) can now be rewritten as 


Therefore 


— = i i = dA , sin 5 
T n T T b " A dt y 


i = tan(e - 6) 


(ri + + i 



^ + csc(e - 5) — 
da 


= tan(e - 5) 


A dt 


y 


(A7) 
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Equation (A7) is the relation between l/R and ^ — that was 

dt 

needed to obtain another relation between the pressure derivatives by 
using equations (A3) and (A4) . If l/R as given by equation (A7) is 
substituted in equation (a4 ) , then 


d j~° s p = -7M 2 tan(e _ 8) 
on 


sin 5 dA , d5 

+ csc(e - 5) — 

y A dt da 


(A8) 


If ^ ^-°S A eliminated between equations (A8) and (A3), then 


dt 


= - 7 M 2 tan(e - 8) 
on 


sin 8 M 2 - 1 8 log p / .v d& 

+ + esc ( € - 8) — 

y 7M 2 dt da 


dt 


A9) 


Solution of equations (A5) and (A9) for the two rates of change of pres- 
sure gives 


d log p _ ylfitan 2 ( 6 - 5) 

dt l - (m 2 - l)t an 2 ( e - 8) 


sin 8 + JL_ cot(e - 5) d log p 
y 7 # sin(e - 5) da 


esc ( e - 8)^- 
' da 


d log p 
dn 


- 7 ^tan( s - 6) 


1 - (M 2 - l)tan 2 ( 6 - 5) 


sin 8 , M 2 - 1 / R \ d log p 

— - — + o- sec ( € - 8) — -- 

y 7M 2 da 


esc ( e - &)— 
da 


By use of the relations 


n = y sec 8 

■p da 

^ de 


and 


d _ d de 
da de da 

_ _1_ _d_ 

Rw de 
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the pressure derivatives can be written 


8 log p 
8t 


7M 2 tan 2 (e - 5) 


1 - (m 2 - l)tan 2 (e - 


_1 n_ cot (e 


5) 


tan 5 + 


5 lQ g P + csc f € _ 

7M 2 Rw sin( 6 - 6) de Rw * 


d& 


(A10) 


d log p _ -7M 2 tan(€ - 5) 

^ log n i _ (m 2 _ i)tan 2 (e - 8) 


tan 5 + 


M 2 - 


jh sec(e - 8) - , j-°g P + — csc(e - 8)|^ 
yW- K w o 6 Ry. d 6 


(All) 


The next step is to find the direction in which dp is zero. This 
direction is that along which the components of the tangential and the 
normal pressure derivatives are equal in magnitude and opposite in sign. 
If 0 p is the angle between this direction and the x-axis, then 

tan ( 0 p - 5) is given by the ratio of p to — log p . Therefore 

' at 8n 


tan(0p - 6) 
tan(e - 6) 


sin(e - 8) sin 8 + — i— cot(e - 8) ^ ^ 

_J_ 7M 2 o € 

-^r sin(e - 8) sin 8 + ^ ~ tan(s - 8) ^ . l°g P + 88^ 

yW- 8 € d e 

(A12) 


The problem that has been solved herein for axially symmetric flow 
has previously been solved by Lin and Rubinov (reference 6) for plane 
flow. Reference 6 also gives the fundamental equations for conservation 
of mass, momentum, and energy in axially symmetric flow. As a check on 
the results obtained in the present paper, equations (A10) and (All), 
the present authors have done the following: By starting with the 

fundamental conservation equations for axially symmetric flow as given 
in reference 6, and by using only algebraic manipulation, as was used 
in reference 6 for the plane case, expressions for the two pressure 
derivatives were obtained. These expressions were identical with those 
derived in the present paper as equations (A10) and (All) . 
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The expressions for pressure derivatives are converted into 
expressions for density and Mach number derivatives by taking into 
account changes in entropy: 


p = e 


AS/c 


v n7 


But 


Therefore 


log P = 7 T^ + 7 log P 
c v 


AS — — (^P — Oy) log 1 i 


log p = -(7 - 1 ) log IT + y log p 


S log p _ 1 S log p 7 - 1 S log IT 
St ~ 7 St 7 St 


Tangent to a streamline, — i is zero. Therefore 

ot 


S log p 



n 


M 2 tan 2 (e - 5) 

(m 2 - l)tan 2 (e - S) 


tan 6 + 


« 


1 

7M“^ 


n COt ( € - S) a log P 
Ry. sin(e - 5) Sc 


U / 

+ CSC ( € 

R w 


- 5) 


S6 

be 


(A13) 
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Normal to a streamline, 

5 log p _ 1 d log p 7-1 5 log IT 

5 log n 7 b log n 7b log n 

= I s lQ g- P + Z--— 1 esc ( c - 5) li£s_n 

7 d log n 7 R w de 

M^tan( e - 5) 

1 - (M 2 - l)tan 2 (e - 6) 


tan 5 + 


M 2 


n 


yM 2 R W 


sec(e - 5) 


d log 


tT esc ( 6 - 6) 
R w 


c)5 

c)e 


+ - — ^ — csc(e - &) 
7 ** 


d log n 
de 


(ai^) 


Equations (A13) and (AlU) give the density derivatives tangent and 
normal to a streamline. Figure 13 shows the variation of the two deri- 
vatives with x/D. 


The angle of the contour of constant density is then given by the equation 
tan ^9p - 


sin(e - 5) sin 5+ cot (e - 8) ^l£S_P + 

7 y-yf oe oe 


tan(e-5) R„. 

— sin(e -5)sin 5 + 

y 


M 2 -! 


7 


tan(e - 8) d lo 8 P , ^5 7-1 1- (M 2 - l)tan 2 (e - 8) d log IT 

M 2 de 7 l^tan(e -8) de 


(A15) 


The variation of 6 with x/D is shown in figure 1^. Furthermore 




7-1 


and therefore 


1 + 


d log M = 


7-1 


yM 2 


M 2 


(d log n - d log p) 


Consequently the two Mach number derivatives are 


d log M 

at 

n 


^1 + - ~ - M^tan 2 ^ - 8) 

1 - (m 2 - l)tan 2 (e - 8) 


1 n cot (e - 8) a log p 
tan 8 + 

yM 2 ^ sin ( g “ 5 ) 


+ csc(e - 8) 


as 

ae 


ro 
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c) log M 
B log n 



7-1 


2 



tan(e - 8) 



l)tan 2 (e - 5) 


tan 8 + 


M 2 - 1 _n_ 
7M 2 ^ 


sec(e - S) 


d log P + 
^6 


n 


CSC ( € 



+ 


2 

7 M 2 


_n 1 5 log TT 

*w sin(e - 8) de 


[\3 

ro 


Then 


tan (e M - 8) 

t an ( £ — 8 ) R^r . / o \ • 

v ' -rf sin(e - 5) sin 8 . _ 

y yMr 


sin(e - 8) sin 8 + cot (e - 8) -■ £ + h 


yW 


+ ^ tan(c-S) — ^ — — - + - — + 


d log p ^8 1 - (m? - l) tan^ ( g - 8) d log IT 

7M^tan(e - 8) 

(A16) 


For use in the final equations (Al4) to (Al6 ) , the value of e for a given value of p/pq 
was found from the equation 


sin^e = 


P/PQ 


7+1)2 (7 - 1] M 2 p 

~r 0 ^ 


Then the corresponding value of x was found from figure 3, y, from figure 2 , and 1^, from 
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figure 5. The value of n is given by n = y sec S. The value of M 
was obtained from the equation 


o o. 2 + (7 - 1) Mo 2 sin 2 e 

Nr - csc^(e - 8) 

27MQ 2 sin 2 € - (7 - l) 


and the value of 8, from the equation 


— Mo 2 

cot 5 = tan e | - 1 1 

lMQ 2 sin 2 e - 1 


(A17) 


From the theory of oblique shock waves, furthermore, 


PO 7 + 1 


27 Mo 2 sin 2 e - 2LlJL 

^ 7 + 1 


(A18) 


n = 


7 + 1 


7+1 

7-1 , 


MQ 2 sin^e 


7 

i7-l 


[7-1 


U + — i MQ 2 sin 2 e 


7Mo 2 sin 2 e - 1 


(A19) 


The derivatives with respect to e that appear in equations (Al4) 
to (Al6) were found by differentiating equations (A17) to (A19) : 


= sin 2 6 
de 


(7 + l)MQ i< 'sin 2 e 


H 1 . h> 8 

+ sec^€ 


^MQ 2 sin 2 e - l) 2 Mo 2 sin 2 e - 1 


d log p 2 7 Mq sin € cos e 


de 


7Mo 2 sin 2 e - 7 g 1 


d log n 27 „ p • * / 1 

— ^ = J-r~T **0 sin e cos e 


MQ 2 sin 2 e 


2 ^ + Mg 2 sin 2 e 


7MQ 2 sin 2 6 - — ^- 1 
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APPENDIX B 

DETERMINATION OF STREAMLINES AND FLOW VARIABLES FROM DENSITY 


The flow field behind the curved shock wave (fig. 2) is rotational. 
The variation in the slope of the shock wave is known (fig. 3)> however, 
as is the density in the flow field (fig. 6). The problem is to deter- 
mine, from these known quantities, the streamlines and the values of 
flow variables other than density. The method that was used was to 
choose an element of arc on the shock wave, to calculate the average 
change in stagnation density across the element of arc, and then to con- 
struct the stream tube step by step in such a way that continuity was 
satisfied. The cross-sectional area of an annular stream tube is 


A = 2rty An 


Continuity requires that the equation 


y An 


p 0 V 0 y 0 (An) 0 
PV 


be satisfied. The velocity V can be expressed in terms of density and 
Mach number by means of the equation 


V = Ma s | 


in which a s y/ a is given by 


*-(!♦ 


7-1 


A 1 ' 2 


in which a s is a constant and in which M is given by 



+ 




1 

7-1 


where p s is a function of the shock- wave angle and free- stream con- 
ditions. The quantity y An can consequently be expressed in terms 
of known quantities in the free stream, the shock-wave angle, and den- 
sity behind the shock wave at An. 


If one side of a stream tube is known, the other side can be 
located. First, the other side is drawn approximately. Then, the 
center line of the stream tube is drawn approximately. With the known 
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value of density and the measured value of y at the intersection of 
the center line and a known contour of constant p/poj the value of An 
is calculated. Then a line perpendicular to the center line is drawn 
through the point. Where this line intersects the known side of the 
stream tube, the value of y is measured. One half of the vertical 
projection of the calculated An is added to this value of y to give 
a closer approximation of y for recalculating An. The new An is 
used to obtain a better approximation for y and so on until the process 
converges. 

The procedure was to obtain first the stream tube adjacent to the 
axis of symmetry and the surface of the sphere. The other stream tubes 
were than located by applying to each in turn the process described. A 
check which showed that any cumulative errors were negligibly small was 
the fact that the method gave deflection angles of the streamlines at 
the shock wave that agreed well with the theoretical values. 
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Figure 1.- Interferogram of flow around sphere at Mach number 1.62. 



Figure 2.- Location of shock wave 
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Figure 3*- Variation of slope of shock wave 
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Figure Variation of second derivative of shock wave. 
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Figure 5.- Variation of radius of curvature along shock wave. 
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Figure 6.- Contours of constant ratio of density to free-stream density 
obtained from interferogram. (From fig. 5(^») of reference 13- ) 
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Figure 7*- Focus of initial tangents of streamlines. 
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Figure 8.- Initial curvature of a streamline 
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Figure 9-- Streamlines and contours of constant Mach number derived from 
continuity and contours of constant density. 
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(a) Contours of constant ratio of density to free-stream density. 


Figure 10.- Results obtained by method of characteristics and by interferometry. 
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(t>) Contours of constant Mach number. 


Figure 10.- Continued. 
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( c ) Streamlines . 


Figure 10.- Concluded 
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Figure 11.- Element of annular stream tube 
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Figure 13.- Variation of density derivatives with x/D. 
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(b) Normal derivative. 
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Figure 14.- Initial angle of density contour as a function of x/D. 
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